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£h . Abstract 

<, 

The Lorcntz contraction of bound states in field theory is often appealed to in quali- 
tative descriptions of high energy particle collisions. Surprisingly, the contraction has 
not been demonstrated explicitly even in simple cases such as the hydrogen atom. It 

■ requires a calculation of wave functions evaluated at equal (ordinary) time for bound 
| states in motion. Such wave functions are not obtained by kinematic boosts from the 

rest frame. Starting from the exact Bethc-Salpcter equation we derive the equal-time 
wave function of a fermion-antifermion bound state in QED, i.e., positronium or the 
hydrogen atom, in any frame to leading order in a. We show explicitly that the bound 
state energy transforms as the fourth component of a vector and that the wave func- 
i tion of the fermion-antifermion Fock state contracts as expected. Transverse photon 

' exchange contributes at leading order to the binding energy of the bound state in mo- 

■ tion. We study the general features of the corresponding fermion-antifermion-photon 
fS |' Fock states, and show that they do not transform by simply contracting. We verify 
^ ■ that the wave function reduces to the light-front one in the infinite momentum frame. 



1 Introduction 

• rH - 

X 

' Bound state wave functions are usually considered only in their center-of-mass frame, 

where rotational symmetry may be fully exploited. In the study of scattering ampli- 
tudes involving several bound states one needs, however, to know the wave functions 
in arbitrary Lorentz frames. Center-of-mass wave functions are commonly defined at 
equal time t = of the constituents, and are non-trivially related to equal time wave 
functions in motion since the relative time is boost dependent £Q. Consequently, the 
Hamiltonian does not commute with the boost generators, causing the boost to be as 
complicated as solving the bound state equation directly in the new frame. 

One manifestation of the non-trivial boost dynamics is the expectation, based on 
classical relativistic physics, that the equal-time wave function of a bound state in 
motion will be contracted in the direction of motion. High energy hadron scattering 
is thus often pictured as collisions between Lorentz-contracted pancakes. This is 
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necessarily a qualitative description since we are far from being able to calculate the 
wave function of a hadron even in the center-of-mass frame. However, it is surprising 
that the equal-time wave functions of much simpler bound states, such as the hydrogen 
atom, have apparently only been considered in the rest frame. 

The hydrogen atom is a non-relativistic state (at leading order in the fine structure 
constant a) and the calculation of its wave function is one of the first exercises in 
courses on quantum mechanics. When the atom is in relativistic motion we must, 
however, make use of the full machinery of relativistic field theory. The excitation 
energies of Fock states with additional particles (electron-positron pairs, or photons) 
may be much less than the energy of the bound state. Hence the contributions of such 
higher Fock states must be carefully considered. As we shall see below, the exchange 
of transverse photons indeed contributes at leading order to the binding of hydrogen 
in motion. This is not unexpected, as a boost of the Coulomb potential leads to a 
transverse electromagnetic field. 

The Lorentz contraction of 3+1 dimensional wave functions in gauge field theories 
has apparently not been demonstrated previously. In P a Lorentz contracting wave 
function of a two body QED bound state is represented as an approximation valid for 
small boosts. The frame dependence of bound state wave functions has been studied 
in various models, see for example, |21E1E1- In El EH Lorentz contraction is obtained 
for a fermion pair interacting via a S potential. There has also been other interesting 
work on the Lorentz covariance of two body equations 0|B]. 

It is reasonable to expect that the equal-time wave function of the hydrogen atom (or 
positronium) can be evaluated analytically in any frame to leading order in a. We 
find that the wave function of the e + e _ Fock state indeed contracts as expected from 
classical relativity. The probability of the e + e~j state is of O (a), which reflects the 
relative scarceness of photon exchange in the weak coupling limit. The photon am- 
plitude does not classically contract, however. More generally, quantum fluctuations 
are unlikely to obey classical transformation laws. 

Rather than trying to boost the well-known rest frame wave function of positronium 
we time order and solve its Bethe-Salpeter equation [5| for an arbitrary momentum 
of the bound state. We generalize our previous 1+1 dimensional calculation |1U| by 
including the transverse photons which contribute in 3+1 dimensions. The equal-time 
formalism is of necessity Lorentz non-covariant - nevertheless we shall see that the 
bound state energy transforms as the fourth component of a Lorentz vector. We study 
the properties of the transverse photon distribution and show that it agrees with the 
light-front result in the limit of an infinitely large bound state momentum. 

2 Wave equation for the hydrogen atom in motion 

We calculate here the equal x° wave function for a (3+1 dimensional) fermion- 
antifermion bound state in the weak coupling limit a <C 1 of QED in any Lorentz 
frame. We will follow closely the procedure of ^U] in 1+1 dimensions: the Bethe- 
Salpeter bound state equation is solved to leading order in a but to all orders in 
|P|/m. Here P is the total momentum of the system and the fermion masses m are 
taken to be equal for notational simplicity. We use Coulomb gauge (k ■ A = 0), but 
also check that the result holds in Feynman gauge. In Coulomb gauge the unphysical 
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Figure 1: The Bethe-Salpeter equation. The blobs represent the wave function \I/p, 
K is the interaction kernel and S is the two-particle propagator. 

photon polarization states are absent and the contribution from physical, transverse 
polarizations is best seen. 

Our starting point, the Bethe-Salpeter equation, is defined as shown in Fig. ^ The 
propagator S is the summed propagator including all radiative corrections. The in- 
teraction kernel K includes all two-particle irreducible interaction diagrams, i.e., di- 
agrams that cannot be split into two interaction graphs just by cutting two fermion 
lines. It is amputated such that it does not include the outgoing or incoming fermion 
propagators which are included in S. The Bethe-Salpeter wave function \I/p is defined 
as the projection of the bound state onto a fermion-antifcrmion state 

*pCp) q /3 = J d*xe**(n\T{#f,(p)1> a (x)}\P\) (1) 

where P is the total momentum, p is the four momentum of the fermion, (f2| is the 
vacuum of QED and A represents all discrete quantum numbers of the bound state 
|PA) (such as spin and orbital angular momentum). 

As in |10| . we work in a time-ordered formalism where the Fock space structure of K 
and \l/p is seen explicitly. The time-ordered rules are obtained by taking a Fourier 
transform over p°. The fermion propagator has forward and backward moving parts 
in (t, p) space 

dp 



S F (t,p) = J ^exp(-itp°) 



■f> — m + ie 

= Q(t) A+(p) cxp{-itE p ) + 8(-t) A - (-p) exp(itE p ) (2) 



where E p = y^p 2 + m 2 and the projection operators A* are defined by 

a^p) ^ ±7 °^ 2 y; p+w . 0) 

Similarly, the photon propagator in (i, p) space reads in Feynman gauge 

(*,p) = -|^[e(t)expH|p|*) + e(-i)exp(i|p|*)] . (4) 

In Coulomb gauge (p • A = 0) we have an instantaneous contribution from the D 00 
component 



Dg»(t,p) = 5(t)±; D%=D% = 
p % p? \ 1 



Dg(*,p) = (^-^j_[e(t)ex P H|p|t) + eM)exp^|p|t)] . (5) 
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Figure 2: Time-ordered single photon exchange diagrams which arise from the inter- 
action kernel K of Fig. ^ ( a ) The instantaneous Coulomb interaction K^. (b) The 
exchange of a transverse photon K^. Time flows to the right. 



In the usual time-ordered perturbation theory one integrates over all time differences 
from zero to infinity. The integrals give energy denominators which are denoted by 
vertical cuts in the Feynman diagrams in the following sections. 

In the center-of-mass frame of positronium the scale of the internal momenta and the 
scale of the binding energy or the potential energy are 

|p| ~ am , AE~V~a 2 m, (6) 

respectively. For a moving system we define the relative momentum by 

q = p P/2. (7) 

We expect the transverse components of q and of the photon momentum k to be the 
same order as in the rest frame ©, whereas the longitudinal components and the 
energy differences will be affected by the contraction, 

q« ~ fcii ~ 'yam , |qj_| ~ |kj_| ~ am , AE ~ 7~ 1 a 2 m (8) 

where 7 = \J~P 2 + (2m) 2 /2m is the boost parameter of the bound state (to lowest 
order in a). 



2.1 Structure of the interaction kernel 

Next we will identify the leading contributions to K using time-ordered diagrams in 
the weak coupling limit. We will see that they arise from the single photon exchange 
part of K which involves only Fock states with one additional photon. The radiative 
corrections to the fermion propagators which are included in S of Fig. prenormalize 
the mass and change the off-shell dependence of the propagator. In a non-relativistic 
bound state the constituents arc nearly on-shcll. Hence to leading order in a these 
effects are accounted for by using the physical mass m in the propagator J5J. 

When iterated the Bethe-Salpeter equation (Fig. gives the wave function as an 
infinite ladder diagram with rungs composed of the kernel K. We will time order the 
ladder and analyze its blocks. We work here in Coulomb gauge, but it is easy to check 
that the results are gauge independent. 



4 



Let us first analyze the single photon exchange diagrams (Fig. [21 . Using time-ordered 
Feynman rules in Coulomb gauge, we have 



K a 

7 



d 3 k i 



E-E p -Ep_ p J (2tt) 3 k 2 

•A+(p)ie7°A + (p - k) <g> A - (P - p + k)ie 7 °A-(P - p) 
d 3 k 1 



E-E p - E P _ p J (2tt)3 2|k| E - £ p _ k - £ P _ P - |k| 1 ] 

SP - J A+(p) Z e7 l A+(p - k) <8> A"(P - p + k)ie 7 J A-(P - p) 

where <g> denotes a direct product between the Dirac spaces of the fermions. For the 
energy denominators arising from the cuts of Fig. 0we have [using JJJ and (JHJ)] 1 

AE F = E- E p - Ep- p ~ AE ~ oPm^ 1 (10) 
AE[ = E-Ep- k -E P -p-\k\=E-Ep-Ep-p + (Ep-Ep^-\k\) 

~ (E p - £ p _ k - |k|) ~ am]" 1 . (11) 

From (|10|) we have for the energy differences AEp ~ aAEj: for the excitation of 
one photon we need an energy O (1/ct) higher than the binding energy in all frames. 
Correspondingly, for the time scales Atj/p = 1/AEjip we have Atj ~ a At p. Hence 
transverse photon exchange is a rare event in the weak coupling limit: The probability 
of finding the bound state in an excited Fock state should be the ratio of the two time 
scales and thus ~ a. This expectation will be confirmed in Sec. |3] 

In Fig. [2 and in the above analysis we assumed that all the fermions move forward in 
time. However, the fermion propagator (J5J also contains a backward moving compo- 
nent. Its inclusion leads to "Z graphs" involving pair production which are suppressed 
in the non-relativistic limit due to large energy denominators. 

Let us pay attention to the coupling structure appearing in . Recalling the defini- 
tion of A* in ©, we use the property {7^,7"} = 2g^ v to write 



A+(p - k) 7 "A + (p) = A+(p - k) 



|^+A-(p)7' i 
Ep 



(12) 



Using A + (p)A (p) = and the results Q, JHJ we can estimate 

\k\/£) = 0(a) 

+ 0(a) (13) 



A+(p-k)A-(p) = 0(\k\/£) = 0(a) 



Ep £ 



where P° ~ £ = ^ (2m) 2 + P 2 . Inserting into we have 

A+(p - k) 7 ^A+(p) = A+(p - k) — + O (a) . (14) 



x We will see that the scaling behavior 7" 1 of AEp is an exact result whereas AEj <x 7 -1 holds 
for fey > 0. For fey < we would have AEj ~ 7am. This reflects the fact that for 7 > 1 backward 
moving photons are suppressed, see Sec. 
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A similar analysis for the antifermion coupling reveals that at leading order we may 
replace 

7 ^A+(p) — — 

A-(p- P )y — -— . (is) 

Using the above estimates we can see that in a general frame both diagrams of Fig. 
13 are of the same order in a. Let us first analyze diagram (a). The first factor in Eq. 
© gives Atp ~ 7 a _2 7n~ 1 . Using (|15l) we have 

A+( P ) 7 WA-(P- P )~ Hp) ~1. (16) 

The Coulomb potential gives using JHJl 

d 3 k 1 k 2 fc 



(2tt) 3 k 2 



a^-/ ~ 7 " 1 a 2 m = U . (17) 



Altogether 

if" - At F ■ V ~ 7a _2 m- 1 • 7 ~ 1 a 2 m = a" (18) 

where we dropped the k dependent projectors A + (p — k) and A~(P — p + k) which 
belong to a different block of K ■ S (see Fig. QJ. For diagram (b) 

~ A+ (P)7 l « V'A-(P - P) ~ (S- -J™)**-^-? (19) 



k 2 j ^ ^ y k 2 ; £ £ r k 2 7 2 

where /3 = |P|/£. The potential contributes 

/" d 3 k 1 i A k^fc,, 

7 2tt) 3 2 k E - E p - £ P _ p+k - k k 



and thus 

7 



^ ~ 7 a~ 2 ™~ 1 • /3 2 7" 2 • 7 a 2 m = /3 2 a° . (21) 



In particular, due to the couplings 1)19(1 . the contribution from transverse photons [Fig. 
[3(b)] is absent at order a in the center-of-mass frame ((3 = 0). This is the expected 
result: for the hydrogen atom at rest transverse photons do not contribute at leading- 
order, but appear as spin dependent interactions at next-to-leading order. The spin 
dependent interactions are hidden in the O (a) terms of <|14|) also when P ^ 0. 

Next we will show that more complicated diagrams can be neglected in all frames. 
A representative set of two photon exchange diagrams is shown in Fig. [21 Diagram 
(a) will be included in our approximation. Diagrams (b) and (d) which include Fock 
states with two photons and diagram (c) will be suppressed. 

Let us study more closely diagrams (a) and (b). Diagram (a) simply consists of two 
separate transverse photon exchanges of Fig. [21 a) , and we have from (|18|) 

~ {K\f ~ /3 4 «° . (22) 
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(a) 



(b) 





(c) 



(d) 



Figure 3: Time-ordered two photon exchange diagrams. Diagrams (a) and (b) arise 
from iterating the covariant Bethe-Salpeter equation (Fig. ^| with one photon ex- 
change kernels. Diagrams (c) and (d) arise from two photon kernels. 

However, in diagram (b) three of the cuts intersect photon lines instead of two. Quan- 
titatively, the only difference to diagram (a) comes from the second cut from the left: 



(a) : 

(b) : 



E — E, 



P k 



£^P-p+k 
1 



7a 2 m 1 



At 1 



E — E,-. — Et>— 



P-p+k+k' 



< 7a m 1 



At j 



(23) 



That is, in diagram (a) the two interactions are separated by the long time scale Atp, 
but in diagram (b) both the interactions must occur within the shorter time scale 
At 1. We have 



K, 



77 



Ati 
At F 



K a 

77 



/3 4 a 



(24) 



and the diagram is thus suppressed. The qualitative picture is that the flight time of a 
photon is short compared to the intervals between the excanges. Thus the probability 
of having two photons at the same instant of time is low. Similar arguments show 
that diagrams (c) and (d) are also suppressed. 

In the case of positronium the kernel K includes annihilation diagrams with, e.g., one 
photon as an intermediate Fock state. Similar arguments as above show that these 
graphs are O (a 2 ) and thus suppressed in all frames. 

We assumed in this section that momenta of order |kj_| ~ oem (and k\\ ~ 7am) 
dominate the integrations, which is true at leading order in a. Transverse photons 
with softer momenta |kj_| ~ a 2 m are suppressed by O (a) due to the smaller phase 
space. However, the flight time of such photons is comparable with the longer time 
scale Atp. This allows an arbitrary number of harder |kj_| ~ am (transverse or 
Coulomb) exchanges while the soft photon is in flight. As the harder interactions are 
O (a ) contributions, the diagrams similar to the one shown in Fig. EUc) but with, 
e.g., several Coulomb exchanges are in fact all 0(a). Such diagrams are known to 
contribute to bound state structure at higher orders in the center-of-mass frame . 
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Figure 4: The time-ordered bound state equation in the ladder approximation. The 
blobs denote the equal-time wave function ipp defined in l|25|) . 



2.2 Lorentz contraction of the e + e wave function 



Now we are ready to write the bound state equation (Fig. ^) in time-ordered form 
with only the leading diagrams included. After the Fourier transform over energy, the 
relevant wave function is the equal-time wave function 



f dp 

ipp(p) aP = / — *p(p)a/3 ■ 

The equation is shown in Fig. 0] The analytic expression reads 

d 3 k 



(25) 



E — E p — Ep-p 



(27T) 



-A+(p) ie 7 V(p - k) ie 7 °A- (P - p) 



E — -Ep_k — Ep-p — |k| E — E p — E] 



P-p+k 



-Ikl 



1 / hi UJ 

2|k| V ~ ~W ) A+(P) ? ^Vp(p - k) iej j A~(P - p) 



(26) 



where ipp is understood as a 4x4 Dirac matrix. In particular, from the property 
A ± A =F = of the projection matrices © it follows that 



A-(pVp(p) = = (^p(p)A+(P - p) 



(27) 



i.e., the wave function l|25(l has only forward moving components in the weak coupling 
limit. 

At leading order in a we may use the replacement (|15l) to eliminate the Dirac structure. 
We have 



<Mp) 



d 3 k 



E-Ep-Ep_pJ (2tt) 3 
1 



^ + ^(P 2 

k 2 8 2 



(P-k) 2 \ 1 



E — En — Ep — 



P-p+k 



|k| E — i? p _k — Ep-p 



5H (28) 



<^p(p-k) 
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E p -E P - p J (2tt) 



,r ' k r y(k)^ P (p-k) (29) 



where £ = \/P 2 + (2m) 2 . The equation thus reduces to a scalar equation for the 
forward moving components of ipp. In comparison with the 1+1 dimensional equation 
of Ref. |1U| . we now have a contribution from transverse photon exchange. 

Let us study the sum E p + -Ep_ p appearing in the denominator of (|29J) . Expanding 
in the relative momentum q of we get 

E p + £ P -p - £ = ^ (qi + 7 ~ 2 <zf ) + O (7 _1 a 4 m) (30) 

where 7 = £/2m and /1 = m/2. Note that the form of ((30(1 is consistent with our 
expectation For the energy denominators in the transverse part of 1)28(1 we get 

E - £ p - k - £ P -p - |k| = (E - E p - E P - P ) + (E p - £ p _ k - |k|) 

P k 



£ 



|k| + 0( 7 - 1 ( 



E - Ep - £ P _ p+k - |k| = -^-|k|+0( 7 - 1 a 2 m) (31) 
and the potential defined in ((29(1 becomes 

—V(k) = 1 1 ^ ( 1 1 1 ^ (32) 

47ra K > k 2 2k 2 ^k 2 + 0k\\ |k| k 2 - /3fc|| |k| J { ' 

- m 1 (33 ) 

k 2 + k2 ( ki+7 _ 2fcf ) 7 2 (ki+7- 2 ^) ( J 

in terms of the parallel and perpendicular components of k and with (3 = |P|/£ . 
In particular, the contribution from the transverse photons [second term in l(32(l ] is 
proportional to 1 and vanishes in the center-of-mass frame. Requiring the potential 
energy 1)33(1 to be commensurate with the relative kinetic energy in 1(30(1 finally verifies 
our expectation (jHJ) and the 7 dependence of the non-leading terms in 130(1 . ((31(1 . 



Let us define the binding energy AM = y/ E 2 — P 2 — 2m which should be independent 
of P. As E — £ = AE ~ 7 _1 a 2 m, we have 

AM = j(E-£)+0 (a 4 m) . (34) 

Inserting (|33J) and (gH) into (21 we get 



Al/-^(qi+7- 2 ?I 



Ana f d a ] 



d 3 k ^p(p-k) 



1 3 k 2 j_ + 7" 2 fc 2 



(35) 



We see that all frame dependence, i.e., the 7 factors, can be removed by rescaling 
fc|| — * 7/c|| and <7|| — * jqu . Then the spectrum is given by the same Schrodinger 
equation as in the center-of-mass frame. The e + e _ wave function exactly Lorcntz 
contracts (or expands in k space) in the direction of motion. Taking into account the 
Dirac structure 1(27(1 . the result can be written 

¥>p(p)a/3 = V ^ x S i, S2 0p(q) (36) 
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where it, v are the usual Dirac spinors and 0p is given by Lorentz contracting the 
usual wave function </>cm of the hydrogen atom at rest: 

<fa>(q) = — ~4>cm (qx, q\\ h) ■ (37) 

The spin wave function \ is normalized as ^2 Sl S9 IXsi,^! 2 = !• At leading order in a 
we were able to replace p and P — p by P/2 in the Dirac structure of l|36|) . 

We finally check that the result (|36J) is also obtained in Feynman gauge with the 
propagator J3J|. In this gauge the photon propagator has no instantaneous part and 
<t2f)l) becomes 

Mp) = E E J_ Ep p J t0 a+ (p) ^ wpCp - k) fer^A- (P - P) 

—1 / z i \ 

2^ U - ^ P -k - ^P-p - |k| + E - E p - £p_ p+k - \k\) ■ { ' 

Using (JTHJl, the D irac structure reduces to a factor —7 2 . Further using CSJ, the 
potential in Feynman gauge becomes 

4^ y ( k ) = ' 2 (k 2 + /3fc|, |k| + k 2 - /3jfc|| |k| ) = _ 7 2 (k^ + 7 - 2 fcj|) ^ 



which is the same result as in Coulomb gauge . The rest of the calculation remains 
unchanged. 



3 The e + e 7 wave function 

In the preceding section we derived the wave function for the e + e~ Fock state, which 
is expected to be the leading component of the bound state in the weak coupling 
limit. Now we will analyze the wave function of the e + e~7 Fock state: we calculate 
the distribution of physical, transverse photons in the bound state. We will see that 
this state occurs with a probability O (a), whereas Fock states with more photons or 
c + c~ pairs contribute terms of O (a 2 ) to the normalization. 

The e + e~7 Fock state contributes as an intermediate state in the derivation of (|35|l . 
see Fig. 0J However, to clarify the derivation of its wavefunction it is best to start 
again from the exact covariant formalism: we will briefly repeat the analysis of the 
previous section but now for the e + e~7 Fock state. The photon distribution will be 
related to the square of the equal-time wave function of the Fock state just as in usual 
non-relativistic quantum mechanics. For an analogous calculation on the light-front 
see [T2G31. 

Let us define the wave function as a coupling to an e + e _ 7 state 

* 7P (p,fc) = J d 4 a;d 4 2/ e+ M -(P- fc ) + ^- fc (O|T{^(0)^( a ;)^(y)} |PA> . (40) 

A covariant equation analogous to the Bethe-Salpeter equation connecting \l/ 7 to the 
wave function ^> of Q may be derived in a way analogous to the derivation of the 
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(a) 




(b) 

Figure 5: Derivation of the Bethe-Salpeter equation for the photon wave function, 
(a) The generating equation for the five-point Green function. G is the four-fermion 
Green function, G 7 has in addition a photon insertion and K 1 is the two-particle 
irreducible kernel, (b) The pole contribution giving the Bethe-Salpeter equation. 
Arrows on the quark lines were dropped for simplicity. 

original Bethe-Salpeter equation (see E])- We only sketch the proof here (see 
Fig. |SJ. The five-point Green function G 1 satisfies the identity shown in Fig. [SJi 
which may be proved diagrammatically. The kernel K 7 is defined to be the sum of 
all two-particle irreducible diagrams, i.e., diagrams which cannot be divided into two 
separate diagrams by cutting two fermion lines. The wave equation is then obtained 
by calculating the residues at the pole caused by the bound state (see Fig. |SJd), i.e., 
at P 2 = M 2 , where P is the total momentum and M is the bound state mass. ^ 7 
is simply given by adding the kernel K 1 to VP . Note that K 1 cannot have a pole at 
P 2 = M 2 : this would lead to G 1 having a double pole. 

Let us move to the time-ordered formalism. The equal-time wave function for the 
single photon Fock state may be defined as 

$£( P ,k) = / * 7P (p,fc) 

= / d 3 xd 3 ye- <x - (p - k) - iy - k <O|#(0)^a(^)^(2/)|PA) 

J x°=y°=0 

= (n|^(0)^ a (p-k)J4"(k)|PA) . (41) 

When time ordered, the interaction K 1 gives rise to several graphs of which some 
are shown in Fig. El Diagram (c) arises in fact as a combination of the lowest order 
kernels K and K 1 , but will be suppressed since the two photons overlap in time. 

Similarly as in Sec. 12.11 we find that diagram (a) together with a similar diagram in 
which the photon is emitted by the antifermion gives the leading contribution in a. 
Thus at leading order in a the wave function l|41(l is related to (p-p of Ij25(l as shown 
in Fig. The analytic expression is in Coulomb gauge 
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(a) 



(b) 



(c) 





(e) 



(f) 



Figure 6: Time-ordered diagrams arising from the interaction kernel Kj. (a) is a 
leading diagram in a. 



P-p 




P-P 




P p+k P p 




p-k 

Figure 7: The calculation of the e + e~7 wave function. 



p-k 



$p(p,k) = 



1 



E - £ p _ k - £p_ p - |k| 2|k| V" k 2 
[^p(p)ie7 J A + (p-k)+A-(P-p)ie7Vp(p-k)] . (42) 



To find the behavior of $ p for the scales of k and q which are relevant in the bound 
state, we may use the previous analysis to simplify the expression. Using (|10[1 and 
l(T5|l we have 



$p(p,k) = 



2k 2 - 2j3fc|||k| V £ 



(^ P (p)-^p(p-k))[l + 0(a)] . (43) 



To give a probabilistic interpretation for the result l|4ri|l we need to define the nor- 
malization of the wave functions. We take the normalization of the bound state to be 
the "non-relativistic" one (without a 28 factor) so that 



(P'A'|PA) = (27r) 3 ^ 3 )(P' - P)S X > 



A • 



Then we have 



(P'A'll, 



(44) 



|PA) (45) 
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where the l's are projection operators on different Fock states 2 , e.g., 
hh = J J^J^J^ 12 c L(Pi)<(P2)4( k )|0)(0|c Sl ( P i)d S2 (p 2 )aA(k) 



Si ,S2 ,A 

V f d ' Pl d3p2 d3k gild 

(2tt)3 (2tt)3 ( 27 r)3 1 1 

• ^( P i)^(p2)^(k)|0)(0|^(p 2 )^( Pl )l M (k) (46) 

in terms of the momentum space field operators defined in (|41|) . Replacing (0| — » 
(Q| and using translation invariance leads to the correct normalization in the weak 
coupling limit 

1= / j^H^^}-J (^(0 2|k|Tr H t ( p ' k ) <i>p ^ p ' k )} + --- 

(47) 

where • • • stands for the contribution from higher Fock states and the trace is over 
the Dirac indices. The probability distribution for the Coulomb gauge wave function 
O (with $° = 0) is thus 



d e V 



(q,k) = ^2|k|Tr{^ t (P/2 + q,k)^(P/2 + q,k)} 

(48) 



d 3 kd 3 q v ' - (2tt) 

a P 2 k 2 ± |0 P (q)-0p(q-k 



4^ 2 |k|3(|k|-/3/fc||) 2 (2^) 3 

expressed in terms of the wave function 0p defined in 1)36(1 . 1)37( 1 . Inserting the dis- 
tribution ((48)1 into the normalization condition 1)47)) and using the estimates © we 
see that the probability for this higher Fock state is O (a). This is in accord with our 
above result that Ati / Atp is O (a) : the fraction of time that a transverse photon 
is being exchanged is of O (a) . It is straightforward to repeat the analysis for even 
higher Fock states and check that their probability is O (a n ) with n > 2. 



3.1 The photon momentum distribution 



The distribution (gHJ may be decomposed into three parts 
a /3 2 k 2 



d 6 V 



d d kd J q 



4vr 2 |k|3(|k|-/3fc„) 
1 



(49) 



(2tt)3 



>(q)r 



;p(q-k)| 2 -2Re{0 P (q) ( /)p(q-k)} 



In the square of the wavefunction ()43)l the first two terms correspond to graphs where 
the photon is emitted and absorbed by the same particle (see Fig. EJ, and the in- 
terference comes from photon exchange. The factor which multiplies the two-particle 
wave functions in 1)49)1 depends only on k and is basically the distribution of photons 



2 The contribution from e+e 7 
state, is not included in 1461 . 



7, which is a higher order correction for a non-relativistic bound 
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emitted by a single fermion. As one might expect, each of the three terms is sepa- 
rately infra-red divergent as |k| — * 0, but the sum is infra-red safe as a consequence 
of the charge neutrality of the system. However, there is an ultra-violet divergence 
as |k| — > oo which stems from the first two terms. For high k the photon wavelength 
is short, and its emission is incoherent. The divergence reflects the distribution of 
photons emitted from a free electron. 

Above we have used an asymmetric coordinate convention because of notational sim- 
plicity. To better understand the distribution, we use in this section a more natural, 
symmetric convention, where the momenta of the fermion, the antifermion and the 
photon are P/2 + q — k/2, P/2 — q — k/2 and k, respectively. That is, we shift 
q — > q + k/2 and thus redefine 



q = p P/2 k/2 



(50) 



in the presence of a photon. We furthermore define the Lorentz contracted momentum 
variables 



k = (kx,fe||) = (kx,&||/7) 
q = (qx,«||) = (q_L ,911/7) • 

In terms of these variables the distribution Ij48(l becomes [remembering (13711 ] 



(51) 



dT 

d 3 kd 3 q 



7/3 a k 



21,2 



4tt 2 



(k 2 L + 7 4|) 3 / 2 (^/kl + 7 2 fef-Ml)' 



9C M 



(q+k/2) 



t>CM 



(q-k/2) 



(2tt) 3 



(52) 



The azimuthally averaged photon distribution, integrated also over the relative mo- 
mentum q of the fermions, reads 



d 2 V 



dk dcos 9 



271 



d^ / d 3 q 



d 6 V 



dk 3 dq 3 

7/3 2 (l — cos 2 1 



4tt 2 
k Jo 



(1 + /3 2 7 2 cos 2 Of/ 2 ( y/1 + /3 2 7 2 cos 2 9 - /3 7 cos 9 



dip 



d 3 q 

(27T) 3 



9CM 



(q + k/2) 



(q-k/2) 



(53) 



where 9 is the angle between P and k, ip is the remaining azimuthal angle and k — |k|. 
Note that the frame dependence only appears in the angular distribution multiplying 
the integrals. An exactly Lorentz contracting distribution would be completely frame 
independent when expressed in terms of the variables H51[). The angular distribution 
is the same as that of a free particle, e — > e + 7. 
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3.2 The photon distribution for the ground state 



The photon distribution l|53|) is valid for any e + e bound state wave function 4>cm- 
Let us study the ground state 



5127T 



1 



(am/2) 2 



(54) 



For the wave function 1|54J) the integral over q in Ij53(l is independent of the angles 9 
and ip. Thus the 9 and k dependence of the distribution completely factorizes, i.e., 



d ' P =^/(cos%(fc) 



dfe dcos 9 2tt 
where the angular dependence is 

7/5 2 (l - cos 2 



/(cos0) 



(1 + /? 2 7 2 cos 2 6») 3 / 2 L/l + /3 2 7 2 cos 2 <9 - /3 7 cos fl) ' 



and 



g(k) ee j 
k 



d 3 q 

(2tt) 3 



^(q + k/2)-^ f (q-k/2) 



(o) 



(55) 



(56) 



(57) 



The dependence of the function / on 7 in (|56J) gives the deviation from exact Lorentz 
contraction. The angular dependence is plotted for various values of 7 in Fig. |S] In 
particular, we have 



/r 2 /(cos#) 

/(cos0) 



(1- cos 2 9) as/?->0 



4 9(cos( 



1 — cos 2 
cos 9 



as P 



1 



(58) 
(59) 



where O is the step function. Thus for large boosts 7> 1 almost all of the photons 
have fe|| > 0. 

Let us study the (frame independent) radial distribution g[k) of (|57|l . For k 3> am, 
the interference term is negligible and we have 



g{k) 



2 

k 



so that the distribution falls as 1/fc for large k. For k <C am we have 



g(k) 



d 3 q 

(2tt) 3 



k-V^(q) 



4/,: 



(am) 2 



(60) 



(61) 



which reflects the decoupling of long wavelength photons from neutral positronium. 
The limit of small and large k behaviors found here is similar for excited states of 
positronium. 

For general values of k the distribution g(k) may be found numerically. It is also 
possible to study the correlation between the magnitudes of q and k by including a q 
dependent cut in the integration of (|57|) . This is illustrated in Fig. Almost all of 



15 



f(cos e) 



/ (y /3"2) 




Figure 8: The angular dependence the contracted and integrated photon distribu- 
tion l|53(l in the positronium ground state. The lines show the angular distribution 
/(cos<9)/( 7 /? 2 ) [defined in ®] for /3 = 0.001, 0.5, 0.9 and 0.999. For (3 = 0.001 (solid 
line) the distribution is close to the symmetric limit (|58|l . For (3 = 0.999 (dotted line) 
the distribution approaches the limit lj59*|) . 




Figure 9: The radial behavior for the integrated photon distribution l|53() . Solid line 
is the function g(k) defined in (|57|) . Dashed and dotted lines are the contributions to 
g(k) from the regions q < 0.7am and q > 2am, respectively. 



l(i 



d A 2 P/dk dcose 



5 - 



A - 



3 - 




Figure 10: The angular dependence of the integrated photon distribution (|62fl for a 
system with velocity [3 — 0.999. Solid, dashed and dotted lines show the angular 
distribution for k = 0.3am, am and 3am, respectively. 



the long wavelength photons k <C am are present when the size of the fermion system 
(given by 1/q) is also large. Conversely, soft photon emission by compact fermion 
systems is suppressed. 

It is also instructive to plot the photon distribution in terms of the usual (uncon- 
tracted) momentum variable k. In terms of k the integrated photon distribution l|53|) 
for the ground state becomes 



d 2 V a I3 2 {1 - cos 2 6) 

dfcdcosfl 2vr (i - /3 C os#) 2 



1 f d 



(2tt) 3 



4 i f (q + k/2)-^°i f (q-k/2)| 2 (62) 



where the integral only depends on k = k^Jl — (3 2 cos 2 9. For (3 — > we have k = k 
and the distribution (f?)2"|l coincides with the k distribution described above. However, 
for large boosts there is a big difference between the distributions. We plot the angular 
distribution for different values of k in Fig. ^|for a system with (3 = 0.999. There is a 
drastic difference when compared to the contracted distribution of Fig. |S] Due to the 
Lorentz contraction effect, the distribution is peaked at forward angles, in particular 
for large k. 

In our calculation (|43|l we assumed that k ~ am. Hence the results do not necessarily 
hold for arbitrarily large values of k. However, it is straightforward to check that the 
1/fc behavior is consistent with the small z behavior ~ 1/z of the Weizsacker- Williams 
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photon distribution function 



a 
2^ 



l + (l-z) 



21 



log- 



(63) 



which gives the probability of finding a physical, transverse photon with a longitudinal 
momentum fraction z in a high energy electron (see, e.g., |15|). 



3.3 High |P| limit 



Finally we check that our result for the photon distribution (I48J1 is compatible with 
the light-front results. It is generally accepted that the usual equal-time picture of 
quantum field theory in the infinite momentum frame coincides with the one calculated 
at equal light-front time x + = t + z = 0. Thus to compare with the light-front results, 
we need to take the limit |P| — * oo. It is easy to check that in this limit our solutions 
for the e + e~ wave function (|3*7jl reduce to the light-front wave functions at leading- 
order in a (see |12l I16| 1. 



»(q) 



'(27T) 



•JLF 



(3/,q±) 



o\ — 
1 IPI 



(64) 



where y is the momentum fraction y = pn/|P| ~ 1/2 + q\\/2"/m. As the two-particle 
Fock state is leading in a also on the light front the normalization can be taken to be 
1 = J dp ± /„ dy\(t) LF (y,pj_)\ 2 . 

The light-front photon distribution is calculated in .12;. As in our calculation the 
probability of the Fock state with one additional photon is O (a). To compare with 
our result we need the relations 



1 



|P|M 



= e(aO- 



l + o 

2|Pb 



m 
IPI 



1 + 



|k|-/3fe|| ~ y ~'{2mx) 2 +k 2 ± 

where x is the momentum fraction x = feii/|P| — ku/2jm. Note that JSJ gives 



y 



Inserting and ||HHJ| into (gHJ we find 



d 3 kd 3 q (q ' k) 



d 6 P 

2 2^ (x, kj_, y, qj_) 

dyd kj^dxd 



(65) 



(66) 



(67) 



kl/x 



\<f>LF{y,q±) - cj>LF{v-x,qx - k±)r 



[(2mx) 2 + k^] 
which agrees with the light-front result |12| 3 . 

3 The result J57J is missing the constraints x < y < 1 which appear at the light front. They are 
unimportant in the weak coupling limit as i ~ a and the light-front wave functions are peaked at 
y~l/2. 
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4 Summary and discussion 

We studied the frame dependence of non-relativistic QED bound states such as 
positronium or the hydrogen atom. Starting from the exact field theoretical Bethe- 
Salpeter equation we evaluated the wave function to leading order in a in all Lorentz 
frames. Using a time-ordered formalism we confirmed the expected result: the e + e _ 
equal-time wave functions exactly Lorentz contract in boosts while the mass spectrum 
is invariant. 

We also solved the leading component of the wave function of the e + e~7 Fock state. 
We saw that this Fock state was a next-to- leading correction with a probability O (a). 
The resulting photon distribution did not contract exactly in boosts, similarly to the 
radiation from a single electron, e — > e + 7. The infinitely boosted limit of the 
distribution was seen to coincide with light-front results. 

The "old-fashioned" time-ordered approach is natural when considering wave func- 
tions evaluated at equal time. The advantages of the time-ordered formalism for 
bound states are well known from studies of positronium in the center-of-mass frame 
- see, e.g., |17j and references therein. In our case the time-ordering was helpful in 
the analysis of the Fock state structure of the bound state and in determining the 
correct order of a for various interaction kernels. 

The covariant (four-dimensional) form of the Bethe-Salpeter equation determines also 
the dependence on the relative time t of the bound state constituents. The time 
dependent wave function \&(f, x) may be derived at weak coupling as in the 1+1 
dimensional case studied in section II. C of EH- The result is 4 



where x\\ = x\\ — ftt is the longitudinal distance between the constituents adjusted for 
the displacement fit = \P\t/E of the center of mass, and (p-p is the equal-time wave 
function ipp of (|36|l Fourier transformed to coordinate space. Combining (|68() with 
our result one may check that the covariant Lorentz transformation formula £Q 
of the e + e~ wave function is satisfied, 



where S(A) is the standard spin transformation matrix of the Dirac equation. 

It would be interesting to extend the analysis presented here to the next-to-leading 
order in a and check if the correction to the wave function contracts classically. 
Likewise one could consider the frame dependence of non-relativistic bound states in 
other theories, e.g., for states bound by scalar exchange. 

In previous work scant attention has been paid to the description of moving bound 
states in field theory, even though much is known about the center-of-mass frame so- 
lutions, in particular for non-relativistic systems. However, the problem is non-trivial 
and certainly worth studying. The usual center-of-mass equal-time wave functions 
are related by an infinite boost to the light-front wave functions which appear in the 

4 The time dependence in the case of non-relativistic center-of-mass motion was given in I5II18I . 




(68) 



* P '(i',x') = 5(A)* P=0 (i,x)5- 1 (A) 



(69) 
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partem model. Hence a better understanding of the behavior of (relativistic) bound 
states under boosts could help to relate the non-relativistic quark model and the 
parton model of hadrons. Evaluating the boosted wave functions of the simpler, non- 
relativistic systems which we studied here is a first step towards understanding the 
boost properties of relativistic systems such as hadrons. 
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